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Recent excperiments (ARPES, Raman) suggest the presence of two distinct energy gaps in high-T c 
superconductors (HTSC), exhibiting different doping dependences. Results of a variational cluster 
approach to the superconducting state of the two-dimensional Hubbard model are presented which 
show that this model qualitatively describes this gap dichotomy: One gap (antinodal) increases 
with less doping, a behavior long considered as reflecting the general gap behavior of the HTSC. On 
the other hand, the near-nodal gap does even slightly decrease with underdoping. An explanation 
of this unexpected behavior is given which emphasizes the crucial role of spin fluctuations in the 
pairing mechanism. 
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The celebrated pscudogap of high-temperature su- 
perconductors (HTSC) is widely believed to be in- 
timately related to the microscopic pairing mecha- 
nism [l|. Early angle-resolved photoemission spec- 
troscopy (ARPES) demonstrated that both the pscudo- 
gap and the superconducting gap below T c are consistent 
with the d x 2_ y 2 symmetry [2J. In addition, both gaps 
were found to have a more or less identical doping de- 
pendence, increasing when the doping was reduced 0, Q • 
This strongly supported a picture, where the pseudo- 
gap is a precursor to the d x 2 _ y 2 -superconducting (SC) 
state, i.e. the pseudogap smoothly evolves into the SC 
gap, when phase coherence of the pairs develops below 
T c [j, @, @| . Accordingly, there seemed to be only one 
energy scale in the system, namely the magnitude of the 
gap in the antinodal region around (7r,0). However, this 
picture has been challenged severely by Raman 0, Q 
and very recent ARPES investigations [E EE [TTJ] : The 
experiments done below T c suggest two distinct energy 
gaps exhibiting different doping dependencies. While the 
gap near the antinodal (jr, 0)-region was still found to 
strongly increase with decreasing doping, the other gap, 
identified from sharp (coherent) peaks near the nodal re- 
gion, has a rather different doping dependence, i.e. docs 
not increase with less doping. These surprising obser- 
vations have been interpreted as a two-gap scenario, in 
which there are two energy scales, one corresponding to 
the near nodal and one to the antinodal region around 
the Fermi surface (FS) . Clearly, these recent findings pose 
also a substantial challenge to the microscopic theory. 

In this letter, we show that the two-dimensional (2D) 
Hubbard model accounts, at least qualitatively, for the 
observed different doping dependencies of the gaps near 
the nodal and near the antinodal regions. We suggest 
(see below) that this "two gaps" scenario can be natu- 
rally explained by the anisotropy of the d-wave coupling 
strength at the Fermi surface, which is related to the 



peaked structure of the spin-mediated pairing interaction 
^Pairing around Qaf = {it, it). These results emphasize 
the importance of spin fluctuations in the pairing mecha- 
nism for HTSC materials El EE EI EE El • In addition, 
they support the premise that the physics of the intrigu- 
ing HTSC materials should be contained in the 2D one- 
band Hubbard model. This premise has recently also 
obtained support from a variety of cluster calculations 
which reproduce the overall ground-state phase diagram 
of the HTSC as well as their single-particle excitation 
spectra [H EE EE HE E0| ■ 

Our main results, obtained by means of the Variational 
Cluster Approach (VCA) [IE EH], are shown in Figs. [2] 
and [3] In Fig. [2j we plot the SC gap near the antin- 
odal and near the nodal point as a function of doping 
for the two-dimensional Hubbard model, whereby the d- 
wave factor has been divided off for a better compari- 
son. The two gaps clearly behave differently as a func- 
tion of doping, in agreement with the experimental find- 
ings [E EE EH- While the antinodal gap increases, the 
nodal gap decreases with less doping. However, an anal- 
ysis of our results in terms of the anomalous self-energy 
indicates that below T c there is only one gap mechanism. 
Both gaps are due to superconductivity, and there is no 
contribution from the normal self-energy, i.e. a remnant 
of the normal-state pseudogap. 

As we will show below, the so-called "two-gap" behav- 
ior can be explained by the doping dependence of the 
spin- fluctuation mediated pairing interaction Vp a i r i ng (q) 
which displays a strongly peaked structure around (ir, n). 
The width a of this structure is given by the inverse 
antifcrromagnetic correlation length and, despite the re- 
duction effect of vertex corrections, a decreases with de- 
creasing doping 0, E. Due to the <i-wave factor, the 
contribution to the superconducting gap Asc(k) has a 
different dependence on the transfer momentum q of the 
effective interaction depending on whether k is close to 
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FIG. 1: (Color online) Left (a): Spectral function of of the 
Hubbard model for U = 8t, t' = -0.3t and x = 0.07 hole- 
doping, calculated with the VCA on a 3 x 3 cluster in the SC 
symmetry-broken phase. The blue arrow denotes the antin- 
odal (AN) region. Right (b): The corresponding Fermi sur- 
face. The antinodal region is marked by AN (blue) , the nodal 
region by N (red). Solid lines indicate the momentum scans, 
along which the gap is determined. Numbers refer to Fig. [2] 



the nodal or antinodal region. This fact, connected with 
the doping dependence of the peak structure of Vp a iring 
around (ix, 7r) accounts for both the qualitatively differ- 
ent doping dependencies of the gaps near (-7T, 0) or near 
as well as for the deviation of the SC gap from the 
nearest-neighbor d-wave form cos k. 
in experiments 0, U , 2(| . 

Our results are obtained on the basis of two numeri- 
cal techniques, which are appropriate for the treatment 
of strongly-correlated materials, namely, the Variational 
Cluster Approach (VCA) 0, 
Carlo (QMC) cluster solution 



cos k y as observed 



23 1 and a quantum Monte 
24, El. The VCA is based 



on the self-energy functional theory, which was proposed 
and applied by Potthoff et al. [23, [H, [2§| . It provides 
a variational scheme to use dynamical information from 
an exactly solvable "reference" system (in the VCA an 
isolated cluster) to go to the infinite-sized lattice fermion 
system at low temperatures and at T = , in particu- 
lar. The ground-state phase diagram of the 2D Hubbard 
model 



U 



E 



(1) 



where Uj denote hopping matrix elements, the den- 
sity at site i with spin "f" and U the local Coulomb 
repulsion was calculated within the VCA by Senechal et 
al. [l8|, and by our group 0, 0, |3(J. For the cluster 
size used in the VCA (up to 10 site clusters solved by a 
Lanczos technique), the T = phase diagram of the Hub- 
bard model in Eq. ([1]), with hopping terms up to third- 
nearest neighbors, correctly reproduces salient features 
of the HTSC, such as the AF and dSC ground states in 
doping ranges, which are qualitatively in agreement with 
electron- and hole-doped cuprates. 

Also the sing le-particle excitations agree with earlier 
QMC data 311 and with the main features observed in 
experiment 32J: clear signatures of the Mott gap, as well 
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FIG. 2: (Color online) (a) "Normalized" SC gap A = 
2Ak/(cosfc a; — cosky) for k near the antinodal (dashed) and 
near the nodal FS point (solid) as a function of doping for the 
Hubbard model (parameters are given in the text), (b) SC 
gap as a function of the d-wave factor (cos k x — cos k y ) on the 
Fermi arc for x — 0.07 hole doping. The numbers correspond 
to the momentum scans marked in Fig. [1] Inset: A(k,u>) in 
a small window around the FS for points to 6. 



as the "coherent" quasiparticle band of width J (= mag- 
netic exchange) and the "incoherent" lower and upper 
Hubbard bands are identified. 

In Fig. [T|b we show a typical Fermi surface for the un- 
derdoped case at x = 0.07 hole doping, obtained from 
the weight of the VCA single-particle spectral function 
at zero energy, A(k, lu = 0). The largest weight is found 
near the nodal direction, marked by N in Fig. [TJ which 
corresponds to the "Fermi arc" . The weak features sug- 
gesting a splitting of the Fermi surface when going to the 
antinodal region are due to the finite cluster used as a 
reference system. Accordingly, we measure the gap in 
the coherent region around the node between line and 
line 6 in Fig. [TJa, and at the antinodal region marked by 
AN in Fig. [TJ 

In Fig. [2k we plot the (normalized) SC gaps A = 
2Ak/ (cos k x — cos k y ) near the nodal and near the antin- 
odal Fermi points as a function of doping x. Like in ex- 
periment, the amplitude of the d-wave gap in the nodal 
region decreases when driving the system from the opti- 
mal doped into the undcrdoped region, while it increases 
in the antinodal region. The different doping depen- 
dences for the nodal and the antinodal regions correspond 
to the two different energy scales observed in photoemis- 
sion experiments. From the figure it is clear that A(k) 
cannot be described by a simple nearest-neighbor d-wave 
form near both the nodal and the antinodal point, es- 
pecially for low dopings. On the other hand, our plot 
of A(k) displayed in Fig. ^Bp shows that in an extended 
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FIG. 3: Spectral function from (ir,0) to at h = 0.07 

doping. Left: A(k, uj) in the SC broken phase. Middle: SC 
phase, but E sc = (see text). Right: Paramagnetic solution. 



region around the nodal point the gap retains its n.n. 
d x 2_ y 2-form, sec also the inset of Fig. [2p. This latter 
finding is again in accordance with experiment (see inset 
of Fig. 3 A in Rcf. 0). 

For the understanding of this behavior it is first of all 
interesting to find out whether the gap originates only 
from one mechanism, related to superconductivity, or 
if there are competing mechanisms producing different 
gaps. In our calculation, we can address this question 
in the following way. In matrix-Nambu formalism the 
Dyson equation reads 



J-T Q(J (Q)-S»°,M -E^Co) 

-E^M w+T a(3 (Q)+E°°*(- 



(2) 

with T Q ^(Q) the hopping matrix, and a, (3 being the 
quantum numbers within the reference system, i.e., the 
cluster sites. After the variational opimization proce- 
dure, we set the SC-related parts of the self-energy S^jg 
to zero and evaluate the corresponding spectral function. 
By comparing the left (full SC solution) and the mid- 
dle panel (S^ = 0) in Fig. [3l one can clearly see that 
the gap around (n,0) is only due to S^g, and there are 
no signatures of a remaining underlying pseudogap pro- 
duced by a mechanism different from superconductivity. 
Fig. [3] shows results for x = 0.07 hole doping, but wc 
have checked that this behavior is the same in the whole 
doping range under investigation. 

Since it is known that above T„ the Hubbard model 
shows a pseudogap behavior 2(| 3l[ , we also compare our 
superconducting results with "normal" -state solutions, 
where we do not allow for U(\) symmetry breaking in 
the variational procedure. The right panel of Fig.[3]shows 
that in this solution there is indeed a pseudogap. Inter- 
estingly, this pseudogap increases with lowering the dop- 
ing and, thus, shows a doping dependence i n q ualitative 
agreement with finite-T QMC simulations [31( and also 
ARPES experiments, identifying the pseudogap above T c 
with the leading edge features near (tt, 0) (JlH, [H|. Our 
calculation is at T = 0. However, taking all the above 
results together, they indicate that the breaking of the 
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FIG. 4: (Color online) The doping dependence of v(q) at 
q = (7T,7r), (f , f ) and (tt, 0) for U = 8t, calculated with QMC 
on an 8 x 8 cluster for inverse temperature f3 = 2.5t~ 1 . Notice 
that the interaction hardly depends on the fermionic momen- 
tum k, so only the dependence on q is shown for k = (— 7r, 0). 



U(l) symmetry destroys the "normal-state" origin of the 
pseudogap, and the gap "below T c " is only due to super- 
conductivity. 

What is the physical reason for the different doping 
dependence of the nodal and antinodal gap ? We argue 
that this behavior is related to the anisotropy of the d- 
wave coupling strength at the Fermi surface in combina- 
tion with the doping dependence of the peaked structure 
of the spin-fluctuation mediated pairing interaction in- 
cluding full vertex corrections. With a previously devel- 
oped formalism [2f| based on QMC, this effective pairing 
strength including the renormalization of the density of 
states can be calculated. Results obtained for the t-t'-U 
(2D) Hubbard model showed that the effective pairing in- 
teraction Vpahing = v(k, q) mediated by AF spin fluctua- 
tions is significantly influenced and reduced due to vertex 
corrections, but it clearly remains peaked around (tt, tt) 



3 l25(. This latter fact is, of course, well known from 
earlier unrcnormalizcd random-phase approximation cal- 
culations. Most importantly, for transfer momenta q 
near (7r,7r), the effective interaction (evaluated at the 
lowest Matsubara frequency) sharply increases with de- 
creasing doping, see Fig.|H In contrast, for transfer mo- 
menta away from (tt, tt), such as q = (tt, 0), v(k,q) only 
shows at best a very weak increase with doping. 

Although HTSC materials are clearly strongly corre- 
lated, our argument can be understood qualitatively by 
means of a simple BCS gap equation for A(kp) at the 
Fermi point : 



A(k F ) 



d 2 k' 
W) 2 



v(k' -k F ) 



A(k') 
E(k<) 



(3) 



Here, E(k') = s/tg, + A(k') 2 , £v being the single- 
particle energy measured from the chemical potential. 
After integrating over momenta perpendicular to the 
Fermi Surface in a thin energy shell of width Q oc J <c t, 
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FIG. 5: Contributions to the gap A(kp) in Eq. [3]as function 
of k', divided by the d-wave factor cos(1cf,x) — cos(kF, y ) for 
normalisation reasons. Typical parameters have been used, 
fi — —1.0t,t' — — 0.3t. Left: Fermi momentum k_F,iv near the 
nodal direction. Arrows indicate transfer vectors with large 
weights. Right: Fermi momentum k_F,Aiv near the antinodal 
direction. The arrow is the AF momentum Qaf- The color 
scale for both plots is given in the right plot. 



one obtains for A<!1: 



A(k F ) Pa - / d% n(k' F ) v(k' F -k F ) A(k' F )log 



A(k' 



(4) 

with n(k' F ) the density of states per unit area of the FS at 
k' F . The integral in Eq. 0] extends over the Fermi surface 
dk' F , the log comes form the integration of the energy 
denominator, and we have exploited the fact that v de- 
pends mainly on the momentum transfer q = k' F — hp. 
For kp near the antinodal point (it, 0), the strongest con- 
tribution (normalized to the value of f(q)) comes from 
k' F near the two antinodal points (0, ir) and (0, — 7r), i. e. 
for q = kj? k' F near (tt, tt). It is thus clear from the 
behavior of v(q) as a function of doping shown in Fig. 0] 
that one expects A(kp ~ (ir,0)) to decrease fast with 
increasing doping. On the other hand, for k F around 
the nodal point (§,§), the contribution for other nodal 
points cannot be strong since it is suppressed by the 
nodes in A(k' F ). As a result, the strength of the SC gap 
is controlled by u(q) at momenta away from q = (ir,ir). 
In order to eludicate this point, we plot in Fig. [5] the 
contributions to the integral Eq.[3]for typical parameters 
fi = — l.Ot, t' = — 0.3t, which show exactly the behavior 
described above. Hence, the behavior of u(q) for q away 
from (-7T, 7r) shown in Fig. 0] provides the physical reason 
for which the SC gap at the nodal point does not increase 
with decreasing doping (33. 34|. 

Summarizing, we have shown that VCA calculations 
within the t-t'-U Hubbard model qualitatively reproduce 
the different doping dependences of the superconducting 
gap recently seen in experiments. When going to lower 
hole dopings, the gap on the Fermi arc, near the nodal 
region, decreases, whereas the gap near the antinode in- 
creases. Moreover we have shown that there is no indica- 
tion in our results for a "two-gap" scenario with distinct 
superconducting and pseudogap below T c . Instead, we 
found that obviously the SC gap in the symmetry bro- 
ken solution absorbs the "normal-state" pseudogap, re- 



sulting in a "one-gap" scenario with only a SC gap below 
T c exhibiting a more complicated structure in momentum 
space. This more complicated structure is shown to be 
naturally explained in terms of a spin-mediated pairing 
mechanism. 
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